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Abstract 

Gauge/string duality is a potentially important framework for addressing the prop- 
erties of the strongly coupled quark gluon plasma produced at RHIC. However, con- 
structing an actual string theory dual to QCD has so far proven elusive. In this paper, 
we take a partial step towards exploring the QCD plasma by investigating the ther- 
modynamics of a non-conformal system, namely the M = 2* theory, which is obtained 
as a mass deformation of the conformal M = 4 gauge theory. We find that at tem- 
peratures of order the mass scale, the thermodynamics of the mass deformed plasma 
is surprisingly close to that of the conformal gauge theory plasma. This suggests that 
many properties of the quark gluon plasma at RHIC may in fact be well described by 
even relatively simple models such as that of the conformal M = 4 plasma. 
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1 Introduction 



The purpose of this paper is to explore properties of strongly coupled nonconformal 
gauge theory plasma within the framework of the gauge theory/string theory correspon- 
dence of Maldacena [1,2]. We are primarily motivated by recent attempts (see [3] and 
references therein) to describe the quark gluon plasma (QGP) produced at RHIC [4-7] 
from the dual holographic perspective. Unfortunately, we do not as yet have a control- 
lable string theory dual to QCD. In fact, most applications of the Maldacena correspon- 
dence to RHIC physics are discussed in the context of conformal Af = 4 supersymmetric 
SU(N) Yang-Mills (SYM) theory in the 't Hooft (planar) limit and for large 't Hooft 
coupling. Remarkably, M = 4 SYM plasma [8-16] as a model for the RHIC QGP 
appears to be reasonably good [17]. We would like to stress, however, that such an 
agreement is rather paradoxical. One requires a strong 't Hooft coupling in order to 
have a controllable string theory dual, and the QGP does indeed appear to be strongly 
coupled. However, it is strongly coupled because it is produced at temperatures of 
order the QCD strong coupling scale, where, naively, conformal invariance is badly 
broken. So, why then does the conformal gauge theory plasma serve as the model for 
RHIC physics? 

As a first step towards answering the latter question, we study the thermodynam- 
ics of the M = 2* (i.e. mass deformed M = 4) SYM plasma over a wide range of 
temperatures and for different mass deformations 1 . The expectation here is that the 
deformation mass scale would provide a model for the QCD strong coupling scale, and 
the thermodynamics of the mass deformed plasma at temperatures of order the mass 
scale would then serve as a more realistic model for the RHIC QGP. 

The paper is organized as follows. In the next section we review the non-extremal 
M = 2* geometry constructed in [21], the holographic renormalization of this the- 
ory [24] and also recall the explicit expressions for the free energy density, energy 
density and the entropy density [24,25]. In section 3 we discuss the map between 
finite temperature M = 2* gauge theory parameters and the parameters of the dual 
non-extremal geometry [26,21]. In section 4 we review the high-temperature thermo- 

1 Although there have been previous attempts to study the thermodynamics of strongly coupled 
nonconformal four dimensional gauge theory plasmas, they have been limited to the high tempera- 
ture regime, where the theory is almost conformal [18-22]. A notable difference is a recent study of 
thermodynamics of strongly coupled Af = 2 gauge theory plasma with massive fundamental hyper- 
multiplets [23]. 
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dynamics of the finite temperature PW flow [21,24,25]. After this, we turn towards 
a numerical investigation of the flow away from the high-temperature regime. We de- 
scribe our numerical procedure in section 5, and present the results of the analysis in 
section 6. Finally, we conclude and also outline open problems in section 7. 



2 Non- extremal J\f = 2* geometry 



The supergravity background dual to finite temperature M — 2* gauge theory [21] is a 
deformation of the original AdSc, x S 5 geometry induced by a pair of scalars a and \ 
of the five-dimensional gauge supergravity. (At zero temperature, such a deformation 
was constructed by Pilch and Warner [27].) According to the general scenario of a 
holographic RG flow, the asymptotic boundary behavior of the supergravity scalars is 
related to the bosonic and fermionic mass parameters of the relevant operators inducing 
the RG flow in the boundary gauge theory. Based on such a relation, and the fact that 
a and \ have conformal dimensions two and one, respectively, we call the supergravity 
scalar a a bosonic deformation, and the supergravity scalar \ a fermionic deformation 
of the D3-brane geometry. 

The action of the five-dimensional gauged supergravity including the scalars a and 
X is given by 



1 -.5 



47TG, 



/ d?V=g [\R - 3(da) 2 - {d X f - V] 
Jm 5 



(2.1) 



where the potential 2 



V = — 
16 



dW 
da 



+ 



dW 
dx 



-W 2 
3 



is a function of a and Xi an d is determined by the superpotential 

1 



W = -e~' la - ^e 4a cosh(2x) . 



In our conventions, the five- dimensional Newton's constant is 

gio . 4t^ 
2 5 voles A^ 2 ' 



(2.2) 



(2.3) 



(2.4) 



2 We set the five-dimensional gauged supergravity coupling to one. This corresponds to setting the 
radius L of the five-dimensional sphere in the undeformed metric to 2. 
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The action (2.1) yields the Einstein equations 



(2.8) 



Rpv = Yld^adyOL + Ad^uX + ^9^V , (2.5) 

as well as the equations for the scalars 

1 &P ldV 

To construct a finite-temperature version of the Pilch- Warner flow, we choose an 
ansatz for the metric respecting rotational but not the Lorentzian invariance 

dsl = -c\{r) dt 2 + cj(r) (dxj + dx\ + dxj) + dr 2 . (2.7) 
With this ansatz, the equations of motion for the background become 

a" + a' (lnc lC 3) -g^ = 0, 
X» + X > f } nc l 4)'-~ = 0, 
c'l + c 1 ! (]n4)' + l Cl V = 0, 
4 + 4 (In Cl c 2 2 )' + ^ C2 V = 0, 

where the prime denotes a derivative with respect to the radial coordinate r. In addi- 
tion, there is a first-order constraint 

(af + I ( X f ~\V- ^(mc 2 )'(mc lC2 )' = . (2.9) 

It was shown in [21] that any solution to (2.8) and (2.9) can be lifted to a full ten- 
dimensional solution of type lib supergravity. This includes the metric, the three- and 
five-form fluxes, the dilaton and the axion. In particular, the ten-dimensional Einstein 
frame metric is given by Eq. (4.12) in [21]. 

For finite temperature flows, we find it convenient to introduce a new radial coor- 
dinate x which spans the range from the horizon to the boundary in finite coordinate 
distance: 

l-x(r) = —, ie[0,l]. (2.10) 

With this new coordinate, the black brane's horizon is at x — 1, while the boundary of 
the asymptotically AdS 5 space-time is at x = 3 . The background equations of motion 



3 This x coordinate is inappropriate for the extremal PW flow, since in that case x is always 
vanishing. In this sense, backgrounds with and without horizons (or cquivalently physics at finite 
temperature and zero temperature) may have distinctly different characteristics. 
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(2.8) become 

4 + 4c 2 (a') 2 



1 



1 



x — 1 



a — 



12Pci(a; - 1) 



x- 1 



x'- 



AVc 2 2 (x-l) 



(x-i){Q( a y + 2( x y) c\ 

- 3c' 2 c 2 - 6(c' 2 ) 2 (x - 1) 

x-1) (6(a') 2 + 2( X ') 2 ) 4 

-34c 2 -6(4) 2 (x-l) 



^ = 0, 



(2.11) 



<9P 

<9x 



0, 



where the prime now denotes a derivative with respect to x. We demand that a physical 
RG flow should correspond to a background geometry with a regular horizon. To ensure 
regularity, it is necessary to impose the following boundary conditions at the horizon: 



x 



a(x), x(x), c 2 (x) 



Sl,S 2 ,S Z f , 



(2.12) 



where the Si are constants. 

The boundary conditions at x = are determined from the requirement that the 
solution should approach the AdS 5 geometry as x — > + : 



x -> 0+ : 



«W> X(x), c 2 (x) 



0, 0, oc x 



-1/4 



(2.13) 



The three supergravity parameters Si uniquely determine a non-singular RG flow in 
the dual gauge theory. As we review in section 3, they are unambiguously related to 
the three physical parameters in the gauge theory: the temperature T, and the bosonic 
and fermionic masses mj and mj of the M = 2* hypermultiplet components. 

A general analytical solution of the system (2.11) with the boundary conditions 
(2.12), (2.13) is unknown. However, it is possible to find an analytical solution in the 
regime of high temperatures 4 (see section 4). Notice that given (2.2) and (2.3), one 
can consistently truncate the supergravity system (2.11) to bosonic deformation only, 
i.e., by setting \ = identically. On the other hand, it is inconsistent (beyond the 
linear approximation) to set bosonic deformation to zero, i.e., to set a = identically 
while keeping the fermionic deformation x ^ 0. 



4 Non-extremal AdS^ geometry is obtained as a trivial solution for vanishing bosonic and fermionic 
deformations: a = \ = identically. 
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In the rest of this section we discuss the asymptotic singularity- free solution of 
(2.11) near the boundary x — > + and near the horizon x — > 1_, constraint by the 
boundary conditions (2.13) and (2.12), respectively. We justify referring to the a and x 
scalars as bosonic and fermionic deformations. We also explain the advantages of using 
the radial coordinate (2.10) in numerical integration. Finally, we recall holographic 
renormalization of the thermal PW flows. 

In what follows we find it convenient to introduce 

c 2 = e A , A = In <5 3 - - ln(2x - x 2 ) + a(x) , 

4 (2.14) 



a 



p = e 

The form of A(x) is chosen to extract the leading asymptotic behavior c 2 ~ x~ l l 4 
from the metric function c 2 (x). In addition, the introduction of p = e a is natural, 
as the bosonic scalar a enters exponentially in the superpotential (2.3). Imposing the 
boundary conditions (2.12) and (2.13), we see that the new warp factor a(x) and scalar 
p{x) satisfy the boundary conditions 



x — > 0+ : a(x) — > , p(x) — > 1 , 
x — > 1_ : a(x) — * In — , /o(x) — > e 1 



(2.15) 



where 5 3 is a new constant. 

2.1 Asymptotics of the thermal PW flow near the boundary 

The most general solution of (2.11) subject to the boundary conditions (2.13) and 
(2.15) may be expanded as a series around x = 0: 

p = 1 + x 1/2 (p 10 + pn In rr) + x (p 2 o + P21 In x + p 2 2 In 2 x) 

(2.16) 




X =Xox 1/4 



Y^Pki Wx H » 

i=l / 

1 + x 1/2 (xio + X11 Inx) + x (X20 + X21 lnx + X22 ln 2 x) 



x - ( 2 - 17 ) 

+ • • • + x fe/2 ( ^ Xfci In 1 x + • • 



,i=i 
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a = x 1 ^ 2 (a w + an In x) + rr (a 2 o + «2i In x + 0:22 ln 2 x) 

*/2 ^ , ^ (2 ' 18) 

+ • • • + x k/2 1 2^ a ^ ln x I + • • • • 

This solution is characterized by five independent parameters 

ln<5 3 ; pn,pio; Xo,Xio \ , (2.19) 



where we have included <5 3 from (2.14). The remaining series coefficients are completely- 
determined from the above parameters. 

For the holographic renormalization (see subsection 2.3) we will need the coefficients 
of the first two subleading terms in (2.16)-(2.18). Explicitly, these are given by 

3 2 

P22 =7^11 , 

26 

P21 =3pi pii - 8p 2 u + yXoPn > ( 2 - 20 ) 

n * 2 3 2 n 26 2 104 2 1 4 

P20 =24p n + -p 10 - 8/01O/0H + yXoPio - — XoPn + ^Xo , 

Xn=^Xo, ( 2 - 21 ) 

_ 9 2 
X22 — -Pn , 

X21 = 2 XoPn + 9pioPn + 7^X0 - yPn » ( 2 - 22 ) 

1 „ 9 2 243 2 51 rt , 13 2 131 4 

X20 =g - 6xoPn + 2P10 + — Pn - yPnPio + 2x Pio + y^XoXio - y^Xo , 



an =0, 



1 2 (2-23) 



«io = - gX 
1 2 

a 22 = - -Pn , 

021 = ~ PioPn - ^Pn - ^Xo , ( 2 - 24 ) 

«20 = - ^PioPn - ^XoXio - yPn + ^Xo - ^P?o ■ 
As indicated above, the non-extremal AdS$ geometry is obtained by setting a = 0, 
X = 0, which in asymptotic expansions (2.20)-(2.24) corresponds to taking pn = 
Pio = Xo = 0. This leads to 

dsl = (2ttT) 2 (2x - x 2 ) ~ 1/2 (-(l-x) 2 dt 2 + dx\ + dx 2 + dx 2 3 ) + -^-^ , (2.25) 
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where the temperature T is given by 

2ttT = 5 3 . 



(2.26) 



Matching with asymptotic extremal AdS$ geometry with a standard Poincare patch 
radial coordinate R (and with our conventional choice for the radius of curvature L — 2) 

R 2 dR 2 
dsj = — (-dt 2 + dx\ + dx\ + dxf) + 4— , (2.27) 

we identify near the boundary 

xocR~\ x-^0 + . (2.28) 

Given (2.28) and the asymptotic expansions (2.16)-(2.18), we identify the conformal 
weight two supergravity scalar p(x) as dual to turning on mass for the bosonic com- 
ponents of the M = 2* hypermultiplet. Parameters p n and p 10 are coefficients of its 
non-normalizable and normalizable modes, correspondingly. Similarly, the conformal 
weight one supergravity scalar x( x ) can be identified as dual to turning on mass for the 
fermionic components of the Af = 2* hypermultiplet. Parameters xo an d Xio are coef- 
ficients of its non-normalizable and normalizable modes, correspondingly. We discuss 
the precise relation of p n and xo to the M = 2* gauge theory bosonic and fermionic 
masses in section 3. Here, we would simply like to emphasize the following well-known 
fact [29]: given {£, Pn, Xo}> the coefficients of the normalizable modes of the super- 
gravity scalars p(x) and x( x )> namely {pio,Xio}j are uniquely fixed by requiring that 
the resulting supergravity RG flow is nonsingular in the bulk. This statement is sim- 
ply the supergravity dual to a gauge theory lore: for a fixed temperature and mass 
parameters, the bosonic and fermionic condensates are determined uniquely 5 . Thus, 
for nonsingular finite temperature Pilch- Warner flows, we must have 

Pio = Pio(4,Pn,Xo) , Xw = Xio(<5 3 ,Pn,Xo) ■ (2.29) 

2.2 Asymptotics of the thermal PW flow near the horizon 

We now consider the behavior of the solution near the horizon, x — 1. Notice that the 
equations (2.11) are invariant under the transformation 

(1 -x) -(1-x). (2.30) 

5 Strictly speaking this is true in the absence of moduli. However, we do not expect moduli at finite 
temperature, and thus unbroken supersymmetry, on the gauge theory side. 
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It is the straightforward to verify that boundary conditions (2.12), i.e. the finiteness of 
{a(x),p(x), x( x )} near the horizon, implies that {a(x),p(x), x( x )} are even functions of 
(1— x) for |1— x\ <C 1. It is this fact that justifies our choice of the radial coordinate x in 
(2.10). Indeed, keeping fixed {£3, pn, xo}, for generic values of {pio, X10}, the functions 
{a(x) , p(x) , xi x )} would diverge at the horizon. Thus, fixing coefficients of the normal- 
izable modes {pio,Xio} from the requirement of nonsingularity of the holographic RG 
flow is equivalent to imposing Neumann boundary conditions on {a(x), p(x), x( x )} a ^ 
the horizon, i.e. as x — > 1_ 

lim p' = lim x' = lim a' = . (2-31) 

x— >1_ X^l- X^l- 

The boundary condition (2.31) would determine bosonic and fermionic condensate 
dependence on the hypermultiplet masses and the temperature (2.29). 



2.3 Holographic thermodynamics of the finite temperature PW flow 



We use the method of holographic renormalization in order to examine the thermody- 
namics of the M = 2* RG flow. Holographic renormalization of this flow was explained 
in [24], where it was investigated using the original radial coordinate r of the metric 
ansatz (2.7). In particular, parameterizing 



Cl 



V) = e A(r) + B(r) ? ^ = ^(r) ? 



(2.32) 



the results of [24] indicate that the entropy density s, the energy density £, and the 
free energy density T are given by 

1 



— — lim e 3A 

t:vJ7Q ^' ^horizon 



1 



lim 

87TG5 r^oo 



J= =£-sT = £ 



3e 4A+B A' + 2e 4A+B {a 1 + a 3 a + a A X 2 + «5 « 2 



+ a & ax 2 + a$ ^— + lne a w x 4 + «nX 4 
me 



1 

8ttG p 



lim 



where 



e = V-9tt 



(2.33) 



(2.34) 
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and the boundary counterterm coefficients ctj take the values 

3 1 1 

on = -, a 2 = - , a 3 = , a A = - , a 5 = 3 , a 6 = , 

3 1 2 1 (2.35) 

«7 = 0, a 8 = --, « 9 = --, «io = -o> «n = ^- 
2 3 3 6 

Using the asymptotic expansions (2.16)-(2.18), and changing the radial coordinate in 
(2.33) following (2.10), we find 



T = -32^ (l + Pn (24 - 96 In 5 3 + 24 In 2) - 24p loPll + 2 X ^io 



+ xS('|-|m2 + |m4 



(2.36) 



Finally, the entropy density is given by 

5 3 e 3a,i 

s = , a h = lim a(x) , (2.38) 

4G5 upl- 
and thus we extract the temperature from (2.37) 

T = — e" 30 " . (2.39) 
2tt v 7 

Of course, the same value of the Hawking temperature (2.39) can be extracted from 
the near-horizon geometry, provided one recalls that sT is a renormalization group flow 
invariant in the supergravity black brane geometries without a chemical potential [30]. 
Indeed, from the relation between components of the Ricci tensor 

R X T = Rt , (2.40) 



we have a constraint 6 



— ) = constant. (2.41) 
C2 ' 



Evaluating the LHS of (2.41) near the horizon, we have 



lim c 4 2 I — ) = 87rG 5 sT . (2.42) 

f' ^horizon \ C2 



3 Eq. (2.41) can also be directly derived from (2.8). 
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Changing to the radial coordinate x in (2.10) and using the asymptotic expansions 
(2.16)-(2.18), the LHS of (2.41) near the boundary takes form 7 

(-£)<4=«- < 2 - 43 > 

Comparing (2.42) with (2.43) leads to the value of sT presented in (2.37). It is pre- 
cisely the holographic RG invariance of the combination sT that guarantees the basic 
thermodynamic relation 

T = E-sT . (2.44) 

Relation (2.44) was verified explicitly for the Af = 2* thermal RG flow in [21]; the 
same relation holds in other thermal gauge/gravity duals, for example in the thermal 
cascading RG flow [22]. 

3 Gauge theory versus supergravity parameters 

The relation between N = 2* gauge theory and the supergravity parameters of the 
(thermal) Pilch- Warner flow was established in [26,21]. For completeness, we review the 
main points here. We begin with the gauge theory, then move to the supersymmetric 
PW flow [26], and finally discuss the gauge/gravity parameter identification at finite 
temperature. 

3.1 N = 2* gauge theory 

In the language of four-dimensional Af = 1 supersymmetry, the mass deformed Af = 4 
SU (N) Yang-Mills theory (Af = 2*) in M 3,1 consists of a vector multiplet V, an adjoint 
chiral superfield $ related by Af = 2 supersymmetry to the gauge field, and two 
additional adjoint chiral multiplets Q and Q which form an Af = 2 hypermultiplet. 
In addition to the usual gauge- invariant kinetic terms for these fields 8 , the theory has 
additional interactions and a hypermultiplet mass term given by the superpotential 

W = ^ Tr([Q, Q]$) + -^—(TtQ 2 + TrQ 2 ) . (3.1) 
9ym 9ym 

7 In principle, one does not have to take a limit in (2.43) since this expression is a constant. 
8 The classical Kahler potential is normalized according to (2/<7y M ) Tr[$<t> + QQ + QQ]. 
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When m — 0, the gauge theory is superconformal with gyM characterizing an exactly 
marginal deformation. The theory has a classical 3(N — 1) complex dimensional mod- 
uli space, which is protected by supersymmetry against (non)-perturbative quantum 
corrections. 

When m ^ 0, the M = 4 supersymmetry is softly broken to M = 2. This mass 
deformation lifts the {Q, Q} hypermultiplet moduli directions, leaving the (N — 1) 
complex dimensional Coulomb branch of the M = 2 SU (N) Yang-Mills theory, param- 
eterized by expectation values of the adjoint scalar 

$ = diag(ai,a 2 , • • • ,a N ) , ^a; = 0, (3.2) 

i 

in the Cartan subalgebra of the gauge group. For generic values of the moduli a*, 
the gauge symmetry is broken to that of the Cartan subalgebra [/(l)^ -1 , up to the 
permutation of individual U(l) factors. Additionally, the superpotential (3.1) induces 
the RG flow of the gauge coupling. While from the gauge theory perspective it is 
straightforward to study this Af = 2* theory at any point on the Coulomb branch [28], 
the PW supergravity flow [27] corresponds to a particular Coulomb branch vacuum. 
More specifically, matching the probe computation in gauge theory and the dual PW 
supergravity flow, it was argued in [26] that the appropriate Coulomb branch vacuum 
corresponds to a linear distribution of the vevs (3.2) as 

di E [-a , do] , al = , (3.3) 

7T 

with (continuous in the large N limit) linear number density 

2 / f a ° 

p(a) — — 5-5 — \/dZ — d 2 , / dap(a) = N . (3.4) 

m 9YM V J-ao 

Unfortunately, the extension of the N = 2* gauge/gravity correspondence of [27,26,31] 
for vacua other than (3.4) is not known. 

In [26, 31] the dynamics of the gauge theory on the D3 brane probe in the PW 
background was studied in detail. It was shown in [26] that the probe has a one complex 
dimensional moduli space, with bulk induced metric precisely equal to the metric on the 
appropriate one complex dimensional submanifold of the SU(N + 1) M = 2* Donagi- 
Witten theory Coulomb branch. This one dimensional submanifold is parameterized 
by the expectation value u of the U(l) complex scalar on the Coulomb branch of 
the theory where SU(N + 1) — > U(l) x SU(N)pw- Here the pw subscript denotes 
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that the SU(N) factor is in the Pilch- Warner vacuum (3.4). Whenever u coincides 
with any of the a« of the PW vacuum, the moduli space metric diverges, signaling the 
appearance of additional massless states. An identical divergence is observed [26,31] 
for the probe D3-brane at the enhancon singularity of the PW background. Away from 
the singularity locus, u = a e [— a ,a ], the gauge theory computation of the probe 
moduli space metric is one-loop exact. This is due to the suppression of instanton 
corrections in the large N limit [26, 32] of M = 2 gauge theories. 

Consider now Af = 2* gauge theory at finite temperature T. Of course, finite 
temperature completely breaks supersymmetry. Thus, we can generalize the thermal 
M = 2* gauge theory by allowing for different (non-equal) masses mj and m,f for the 
bosonic and fermionic components of the M = 2* hypermultiplet {Q, Q} correspond- 
ingly. It is only when m& — m.f — m, and at zero temperature T = 0, that we have 
M = 2 supersymmetry. Since turning on mass terms for the bosonic or fermionic 
components of the hypermultiplet corresponds to deforming the M = 4 supersym- 
metric Yang-Mills Lagrangian by relevant operators of different classical dimension (a 
dimension two operator for bosonic and a dimension three operator for fermionic mass 
terms), such deformations will be encoded in different supergravity modes. As indi- 
cated in section 2, turning on bosonic/fermionic masses corresponds to turning on the 
five-dimensional supergravity scalars a/x correspondingly. 

We would like to conclude this section with a simple observation. Turning on 
bosonic/fermionic masses for the hypermultiplet components sets a strong coupling 
scale A oc max{mb,m/}. In this case, we expect to find two qualitatively different 
phases of this gauge theory, depending on whether T > A or T C A. In the former 
case we expect the thermodynamics to be qualitatively (and quantitatively in the limit 
T/A — > oo) similar to that of the M = 4 gauge theory plasma [33]. On the other 
hand, when T ~ A, and with mj = 0, we actually expect an instability in the system. 
Indeed 9 , turning on only the supergravity scalar a, i.e., setting m& ^ and rrif = 0, 
corresponds to giving positive mass-squared to four out of six M = 4 scalars (these 
are the bosonic components of the M = 2 hypermultiplet). The remaining two M = 4 
scalars at the same time obtain a negative mass-squared [27] — they are tachyons at 
zero temperature. At high enough temperature the thermal corrections would stabilize 
these tachyons. However, as we lower the temperature, we expect the re-emergence 
of these tachyons. As argued in [34], dynamical instabilities in thermal systems are 

9 We would like to thank Ofer Aharony for clarifying this. 
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reflected in thermodynamic instabilities. Furthermore, it was argued in general (and 
demonstrated explicitly with a concrete example) [35] that thermodynamic instabilities 
are reflected to developing (? s < 0, where c s is the speed of sound waves in the thermal 
gauge theory plasma. We discuss the supergravity realization of these phenomena in 
section 6. 



3.2 Supersymmetric PW flow 

It is clear that the following conditions must hold in order to preserve J\f = 2 super- 
symmetry for the mass deformation of J\f = 4 SYM: 

■ the temperature must be zero: T = 0; 

■ the masses for the bosonic and fermionic components of the Af = 2 hypermultiplet 
{Q, Q} must be the same: m h = mj = m. 

The former condition corresponds to a restriction 

ci(r) = c 2 (r), (3.5) 

ensuring Lorentz invariance of the metric. In this case, the supergravity RG flow 
cannot be parameterized by the radial coordinate x introduced in (2.10). Instead, the 
supersymmetric flow is easiest to parameterize in terms of the fermionic scalar \ 
[0, +oo), where x = corresponds to the asymptotic AdS§ boundary, and \ ~^ +°° to 
the enhancon location in the bulk [26,31]. The supersymmetric Pilch- Warner solution 
is then given by [27] 

e A - ^ 



sinh(2x) ' ^ ^ 

p 6 = cosh(2 X ) + sinh 2 (2 X )ln^4, 

cosh(x) 

where the single integration constant k is related to the hypermultiplet mass m ac- 
cording to [26] 

k = mL = 2m . (3.7) 

In order to identify the thermal RG parameters {pn, Xo} with the gauge theory masses, 
we need the boundary asymptotics of the PW flow. This was computed in [21]. Intro- 
ducing 

x = e- r/2 , (3.8) 
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we have 

X = kx 1 + k 2 x 2 (| + | \n(kx)) + k A x A (-£ + f \n{kx) + f \n 2 (kx)) 
+ 0(k 6 x 6 ln 3 (kx)) 

p = 1 + A; 2 £ 2 (I + | \n(kx)) + k A x 4 (i + 2 ln(fcr) + § ln 2 (A;£)) + (A; 6 £ 6 ln 3 (A;£)) , 
A = - ln(2x) - \k 2 x 2 - k 4 x A (| + f ln(ifex) + § ln 2 (A;£)) + C (k 6 x e ln 3 {kx)) . 

(3.9) 

3.3 Thermal PW flow 

We now consider the thermal PW flow with m& = to/ = to. Here we would like to 
identify the parameters pn and Xo f° r sucn a fl° w i n terms of the PW quantity k 
(and thus, following (3.7), the hypermultiplet mass to) by comparing the boundary 
asymptotics (2.16)-(2.18) with (3.9). In order to do so, we first establish the relation 
between the radial coordinate x in (2.16)-(2.18) and the radial coordinate x in (3.9). 
Matching the asymptotic warp factors e A for the zero temperature PW flow and the 
thermal PW flow as x — > 0+, we find 

x~85 4 x 4 , x^0 + . (3.10) 

Next, using (3.10) and matching the leading nontrivial asymptotics [order x l / 2 hix for 
p(x) and order x l / A for x( x )} f° r both the thermal PW solution and the supersymmetric 
PW flow, we find 

2 1/2 k 2 V2 6n /m\2 2 1/4 k 1 , a fm\ 

= = e (r) • *> = tit = ^ e • (3 ' u) 

where we used (2.39). We emphasize that the above matching procedure is well- 
defined — it was shown in [21] that finite temperature effects modify the asymptotic 
supersymmetric PW geometry at order 0(x) ~ 0(x A ), which is subleading compared 
to the order at which the matching (3.11) is done. 

Motivated by (3.11), we now propose that in the general case, i.e. for m b ^ rrif, we 
may independently extract bosonic and fermionic masses according to 

on - — e" 6 ^ (^Y Yn - — e-^ (3 12) 

In section 4 we recall that such an identification leads to a consistent thermodynamics 
of the thermal PW flows [24, 25] at high temperatures. Furthermore, we later show in 
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section 6 that this identification in fact leads to a consistent thermodynamics at any 
temperature. 

In the rest of this subsection we derive some useful expressions for the free energy 
and the speed of sound waves. Specifically, we would like to express the free energy 
and the speed of sound in terms of the coefficients of the non-normalizable modes 
of the supergravity scalars p(x),x(x). We also present the differential constraint on 
coefficients of the normalizable modes of these supergravity scalars following from the 
first law of thermodynamics. 

Since we will constantly refer to the parameters pu and Xo related to the bosonic and 
fermionic masses, we distinguish them from derived quantities (such as the subleading 
coefficients) by introducing the notation 

f = Pn, V = Xo- (3-13) 

Notice that because of (3.12), £ and rj are not independent, but are related through 
the ratio of bosonic and fermionic masses according to 

V 2 = Q^Z- (3.14) 

m b 

The nonsingularity of the thermal PW flow at the horizon (2.31) will determine the 
coefficients {pio, Xio} of the normalizable modes of the supergravity scalars {p(x), x( x )} 
and also the value a h of the warp factor a(x) at the horizon in terms of (3.13) and 
(3.14). When m b ^ 0, we take 

Pio = Pio(0» Xio = Xio(0; ah = ah(0- ( 3 -15) 
On the other hand, for m b = 0, we have instead 

pw = pio(v), Xw = Xio(v), a h = a h (rj). (3.16) 

Notice that a given pair {£, r]} can be unambiguously related to bosonic and fermionic 
masses, measured with respect to temperature. Indeed, following (3.12) 

(^) 2 = 12 V2 tt 2 e ** £ , (^) = 2 3 / 4 tt e ** rj . (3.17) 

Solving for 5^ from (2.39) and expressing G5 in gauge theory variables (2.4), we find 
for the free energy 

T = - Wt 4 e ^(0 (l + 24£ 2 ln(£ 2 ) (l " ^) + ^ ^ Xio(0 - ^ Pio(o) 
+ ^0, 

(3.18) 
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assuming m b ^ 0, and 

? = - \n 2 N 2 T* e 12a ^ (l - \rf ln^ 4 ) + 2r? 2 Xl0 (r ] ^j + T , (3.19) 

if £ = 0. In the above expressions, JF does not depend on temperature. Such a constant 
arises because holographic renormalization generically 10 leads to In T dependence in the 
free energy, which requires the introduction of an arbitrary regularization scale /x 

InT — > In — . (3.20) 

Changing this scale modifies JF , but otherwise has no effect on the thermodynamics. 
This phenomena has been encountered previously; it was discussed in perturbative 
gauge theory thermodynamics in [36], and in the context of gauge/string theory cor- 
respondence in [37] 11 . Notice that in (3.18) we wrote ln(£ 2 ) instead of 21n(£). The 
reason is that the former expression allows us to study thermodynamics for £ < 0, 
which according to (3.12) we interpret as m 2 < 0. While introducing tachyonic masses 
for gauge theory scalars at zero temperature leads to instability, this is not the case 
at finite temperature. Here, the effective mass squared receives thermal corrections of 
order T 2 , which might cure the zero temperature instability. 

Besides the free energy, we can also evaluate the entropy density and the energy 
density 

s = ^ti 2 N 2 T 3 e 12ah , £ = F + sT. (3.21) 

In addition, from the defining equations (3.17), we can evaluate the quantities d^/dT 
and drj/dT: 



dT 6V2n 2 T* \ <% J 



dT 

dr l _ m f „-3a h („) f „ da h^) 



(3.22) 

Using (3.18) or (3.19), (3.21) and (3.22) the first law of thermodynamics, 

df=-sdT, (3.23) 



10 Notice that this is not the case when m\> = mf. 

11 It appears even in the high temperature thermodynamics of the thermal PW flow [25] (see also 
section 4). 
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reduces to a differential constraint on the quantities 12 

{£,Xio(0>Pio(0> a /i(0} form 6 ^0; 
{v,Xw(v),Pw(v),ah(v)} form 6 = 0. 
In particular, the first law in differential form reads 



(3.24) 



=4£ (l-5) + 2 (p 10 (0 - C-^ ) - ^ ( X 10 (0 - ^ ) + 



(3.25) 



for m b ^ 0, and 



= - V +2 V xM - rj - 6— , (3.26) 



for nib = 0. Using either (3.25) or (3.26), we can evaluate the speed of sound entirely 
as a function of £ or rj. Straightforward algebraic manipulations result in 



, dP dT 1 + E 



° s ~ dS ~ d£ ~ 3 - E ' (3 ' 27) 



where 

,4 



E = E(0 =24£ 2 - 1 + 12$ - Pio(0 



m 6 V V s 

2 



+ «e^ (xi»(f)-{%^), (3.28) 

v _ v./ x 2 4 . 2 / \ 1 3 d Xw(v) n 
£ = £(77) =-77 + 77 - g ^ dr] , rn b = 0, 

and where P = —T is the pressure. 

In general, one may explore the PW thermodynamics in the entire phase plane of 
m b -mf. However, here we focus on two special cases: 
■ when nif — (we call this the bosonic case) 13 we have 

n 2 N 2 T 4 



bosonic 



e 12ah{& (l + 24£ 2 ln(£ 2 ) - 24£ p 10 (o) + ^0 , (3.29) 



Ef_(0 = -24£ 2 - 12£ Pl0 (O + 12£ 2 , (3.30) 



12 Notice that even though, for £ = 0, the free energy density (3.19) does not depend on p w , one 
generically expects that pw = pio(j]). 

13 Recall that this corresponds to a consistent truncation of the supcrgravity 
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and the first law of thermodynamics takes form 

= ^ + 2(p 1 o«)-«^)+^f ; (3.31) 

■ when rrif = rrib = m (we call this the supersymmetric case), or correspondingly from 
(3.14) 

V 2 = 6Z, (3.32) 

we have 

Tsusy = g e 12a ^ f 1 - 24£ pi (O + 12f Xio(0 J + , (3.33) 

WO = 12C - PIO(O) + 6^ (x 10 (0 - , (3.34) 

and the first law of thermodynamics takes form 

= 2 (p 10 ({) - {*^) - ( Xl „ K ) - f*^0) + ^> . (3.35) 

Notice that in the supersymmetric case there is no ln(£ 2 ) dependence in the holographic 
free energy (3.18). As a result one can unambiguously determine JF directly from 
(2.33) 14 . We find 

8 S 288tt 2 v ; 

4 Thermodynamics of the high temperature Pilch- Warner flow 

We now recall that the differential equations (2.11) describing the finite temperature 
PW renormalization group flow admit a perturbative analytical solution at high tem- 
perature [21,24,25]. In this regime, the appropriate expansion parameters are 

5iOc(^) 2 «l, 5 2 oc^«l. (4.1) 

In this case, the metric function A(x) is given to first nontrivial order in Si, 82 by 

A(x) = \n5 3 - i In (2x - x 2 ) + 5f A 1 (x) + b\ A 2 (x) , 

a(x) —81 ati(x) , ( 4 - 2 ) 
X(x) =8 2 X2(x) , 



14 The latter holographic prescription was formulated in such a way that the holographic energy of 
the supersymmetric PW flow vanishes [24] . 
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where 



ai = (2x-x 2 ) 1/2 2Fl {\,\-l-{l-xf) , 
X2 = (2x-x 2 )^ 2 F 1 (H;l;(l-x) 2 ) , 



A x =4 
A 2 4 



1 (z - l)dz 
x (2z-z*)* 

4 f 1 (z-l)dz_ 

2Y2 



( 7l " I dy 



daA 2 (2y - y 
dy J y-l 



2\2 



3 J x (2z - z 



72 



dy 



2\2 



d X 2Y (2y-y 2 ) 
dy J y-l 



(4.3) 
(4.4) 

(4.5) 



The constants 7$ are fine-tuned to satisfy the boundary conditions, and are given by 

8 - 7t 2 8 - 3tt 



7i = 



72 



2tt 2 ' " 8vr 
Comparing (4.3), (4.4) and (2.16), (2.17) we find the relations 

__V2 

71 

2^0F 



(4.6) 



P11 



3\ 2 



r(!) 



*2, 



3V2 ln2 . 
P10 = z ^1 > 



7T 



X10 = 



2 l/2 r( 3)4 



(4.7) 



7T 



Also notice that 



2ttT =5 3 1 + 5 2 2 - 



<ir 2 



1 d 2 A 2 



<ix 2 



x=l 



7T Z 671 



(4.8) 



The parameters {^i,^} can be related to the {m b /T,mf /T} parameters of the dual 
gauge theory via 



Si 



1 



m b \ 2 
24tt V T J ' 



[r(»]' 



(4.9) 



2tt 3 / 2 T ' 

To leading order in mass deformation, the free energy (3.18) takes form [24,25] 



—tt 2 N 2 T 4 
8 



192 



1 - — ln(TrT) 5( - - 51 

7T Z TT 



(4.10) 



One can further evaluate the entropy density of the non-extremal PW geometry, 

s = ViV 2 T 3 ( 1 - 5\ - 1 5 2 ) , (4.11) 
and verify that the first law of thermodynamics (3.23) is satisfied. 
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5 Numerical procedure and results 

We are mainly interested in the thermodynamics of the PW flow away from the high 
temperature limit. In this case, since analytic solutions to the equations of motion 
(2.11) are unavailable, we resort to numerical integration. In this section we describe 
the numerical procedure and some of the consistency checks we have performed on the 
numerical results. Extracting the actual thermodynamics from these solutions will be 
taken up in the following section. 

The numerical procedure is conceptually the same for both the bosonic and su- 
persymmetric mass deformations of the non-extremal AdS§ geometry. In practice, 
however, it is much easier to implement for the bosonic deformation 

to/ = , TOf, 7^ , (5.1) 

and so we will describe the details of this first. We then outline the modifications 
required for handling the supersymmetric mass deformation 

rrif = m b = to , (5.2) 

and comment on thermal flows corresponding to generic mass deformations. Finally, 
we discuss the validity of the supergravity approximation. 



5.1 Bosonic deformation 

Simplicity of treating the bosonic mass deformation stems from the fact that thermal 
PW flows allow for a consistent truncation with the supergravity scalar x( x ) identically 
set to zero. The relevant equations can thus be obtained from (2.11) with 

x(x) = , a = \np, In c 2 (x) = A(x) = In <5 3 — - \n(2x — x 2 ) + a(x) . (5.3) 

The boundary conditions (which guarantee that the flow is singularity-free) are given 
by (2.15) and (2.31) 

x — > 0+ a(x) — > , p(x) — > 1 ; 
x^l_ a'(a;)->0, p'(x)^0. 

Asymptotic expansions for {p(x), a(x)} as x — > + are given by (2.16), (2.18). The first 
couple terms are presented in (2.20)-(2.24), where one has to set Xo = corresponding 
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to xfa) = 0) see (2-17). Notice that numerical integration of pfa) and afa) from 
x = + is uniquely determined specifying the non-normalizable pu and normalizable 
Piq mode coefficients. However, a generic choice of a pair {pn,pio} would produce a 
solution that does not satisfy the horizon boundary condition (the x — > 1- boundary 
condition in (5.4)). Solving Sturm-Liouville problem with boundary conditions (5.4) 
would determine 

Pio = Pio(pn) Pio = Pio(0 , (5.5) 

where on the RHS we used conventions of subsection 3.3. For a flow satisfying (5.5) 
we can extract 

a h {£) = lim a(x) . (5.6) 

X— »1- 

To summarize, each non-singular bosonic deformation flow would generate a triplet of 
numbers 

{£,Pio(£WO}, (5-7) 

which can be used to evaluate the free energy J-'bosonic (3.29), the speed of sound c 2 s bosonic 
(3.27) and (3.30), and verify (numerically) the first law of thermodynamics (3.31). 

We use the following algorithm to generate triplets (5.7). First, we choose a mass 
deformation parameter pn = £ which we keep fixed during the iteration procedure. 
The iteration starts by choosing a trial value of pio- Given {pu, pw} we integrate 
numerically 15 (2.11) (with (5.3)) from x ini tiai < 1 to x final (such that (1 -x final ) < 1) 
using the power series (2.16), (2.17) to specify 16 

P fainitial) i P fa initial) i & fa initial) 5 ^ fainitial) ■ (5-8) 

For a generic value of p w , 

lim p fa) — > ±00 , lim a fa) — > — oo . (5.9) 

We find that for a given value pu = £ there is (generically) a unique p 10 = Pio(0> such 
that pfa) and afa) are finite as we go to the horizon, x — > 1_, i.e., 

J™ ~^faiP u = ^P 10 = Pio(O) = 0, ( 5 - 10 ) 



15 We used Mathematica for numerical integration as well as our own C code based on fifth-order 
Runge-Kutta. Both procedures (up to controllable numerical errors) produced equivalent results. 

16 In practice we truncated the power series to k = 3 terms (inclusive), and used x initial = 10~ 9 , 
Xfinal = 1 - 10~ 6 . 
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Pio(0 



-0.2 -0.15 




Pcrit 



ah{i) 



-0.2 -0.15 -0.1 




crit 



Figure 1: Altogether we obtained 3227 triplets {£, pio(0' a h(0}- The left plot rep- 
resents piolO; an d the right plot represents a/i(£). The critical value of £ is at 
p crit w 0.03582(2). Red segments of the plots correspond to £ G [0,p cr it) such that 
Pio(0) = a/i(0) = 0. Blue segments of the plots do not pass through the origin and 
correspond to < £ < p cr u. Black segments of the plots correspond to £ < 0. 



while 



£,Pio^Pio(0) = ±°o- (5-11) 

Now, the iterative algorithm is clear: we use Newton's shooting method of determining 
pio(0 for a fixed p u = £. 

Results of the numerical analysis are presented in Figs. 1-3. From Fig. 1 notice 
that for pu > p cr it ~ 0.03582 there is no solution to Sturm-Liouville problem with 
boundary condition (5.4). Using (3.17) we find that corresponding to p cr u, 



T 



2.29(9). 



(5.12) 



crit 



We find that for each value of £ G (0, p cr it) there are two sets {pio(0> a h(0} which are 
represented by red/blue segments of the plots in Fig. 1. For each £ < we find a single 
set {pio(0> °/i(0} represented by black segments of the plots in Fig. 1. In section 6 
we show that red/blue values of the bosonic condensate pio(0 correspond to different 
phases of the M = 2* plasma. We also show that there is £ = p* < p cr u such that for all 
£ > p*, M = 2* plasma becomes unstable both thermodynamically and dynamically, 
in agreement with general arguments of [35]. 

Figs. 2-3 represent p'i (C,), c4(£), and a numerical verification of the first law of 
thermodynamics (3.31). Since our numerical data is rather densely spaced 17 in £, we 
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Mostly we have A£ = 10" 
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use 

/(0 = " 5r = , (5-13) 

as a numerical approximation for a derivative of a function /(£)• Notice that 

hm — — — = =foo , hm — — — = =foo , (5-14) 

with minus/plus signs for red/blue segments of the curves in Fig. 2/Fig. 3 correspond- 
ingly. Also, from Fig. 3 

hm — — — = +oo , hm — — — = +oo . 5.15 
i^+o d£ e^+o d£ 



We find that 



MO 



« -2.10(6), (5.16) 

5=+2xl0- 4 

which is rather close to a predicted value from the high temperature analysis of section 
4 (see (4.7)) 

high temperature 

p' 10 (0 = -3 In 2 « -2.07944. (5.17) 

A better agreement is achieved if we fit the first 200 points of the red segment of pio(£) 
on the left plot in Fig. 1 with a 10th order polynomial in £. We find in this case 



p'io(0 



fit 

^-2.07972. (5.18) 

5=o 



The green points on the right plots in Figs. 2-3 demonstrate the cancellation of a' h (£) 
with a combination 

as required by the first law of thermodynamics, see (3.31). We find that such a cancel- 
lation is achieved with an accuracy of better than 10~ 4 . Clearly, this provides a rather 
impressive check on both the validity of the holographic renormalization of the thermal 
PW flows explained in [24], and the numerical procedures developed in this paper. 

5.2 Supersymmetric deformation 

Albeit more complicated, the Sturm-Liouville problem for the supersymmetric mass 
deformation is conceptually similar to the one for the bosonic mass deformation studied 
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P'io(0 a' h (0 




Figure 2: The left plot represents p' 10 (0 for both the red and the black segments of 
the left plot in Fig. 1. The right plot represents a' h (C,) for both the red and the black 
segments of the right plot in Fig. 1. The green points on the right plot represent 
numerical verification of the first law of the thermodynamics, where a' h (£) must be 
canceled with an appropriate combination of {£, pw(£,), P'(Oh see (3-31). 

in previous subsection. Here we have to use all equations in (2.11). As for the bosonic 
deformation, we rewrite these equations using {p(x),a(x)} such that 

a = lnp, \nc2(x) = A(x) = In 53 — — ln(2x — x 2 ) + a(x) . (5.19) 

The boundary conditions (which guarantee that the flow is singularity-free) are given 
by (2.15) and (2.31) 

x — > 0+ a(x) — > , p(x) — > 1 , y(x) — > ; 

k j pv j (5.20) 

x -> 1_ a'(a;)->0, //(z)->0, x'(^)^0. 

For the bosonic mass deformation we have a one- dimensional Sturm-Liouville problem: 
given pn, the boundary condition at the horizon determines p w ; the pair {p n = 
C)Pio(O) further determines the value of the a(x) at the horizon, For the 

super symmetric mass deformation the Sturm-Liouville problem is two-dimensional: 
given pn = £, the boundary condition (5.19) would determine 

Pio=Pio(0» Xio = Xio(0- ( 5 - 21 ) 
For a flow satisfying (5.21) we can extract 

a fc (f) = lim o(z) . (5.22) 
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Figure 3: The left plot represents p' 10 (0 for the blue segment of the left plot in Fig. 1. 
The right plot represents a' h (£) for the blue segment of the right plot in Fig. 1. The 
green points on the right plot represent numerical verification of the first law of the 
thermodynamics, where a' h (£) must be canceled with an appropriate combination of 
{£,/>io(£), //(£)}, see (3.31). 

Thus, each non-singular supersymmetric mass deformation flow would generate a set 
of four numbers 

{£,Pio(0,Xio(0,^(0}, (5.23) 

which can be used to evaluate the free energy T susy (3.33), the speed of sound <? s susy 
(3.27) and (3.34), and verify (numerically) the first law of thermodynamics (3.35). 

Since the Sturm-Liouville problem is technically more complex here, we restrict our 
discussion to £ > 0, or following (3.12) m 2 > 0. Results of the numerical analysis are 
presented in Figs. (4)- (5). Fitting the first 200 points of pio(£) and Xio(0 with a 10th 
order polynomial in £ we find 



P'io(0 



fit 



5=0 



-2.069(2) , 



Xio 



fit 



5=0 



-0.3219(7) 



(5.24) 



which should be compared with the high temperature predictions 18 (4.7) 

high temperature 



P'io(0 



Xio(0 



high temperature 
?=0 



-3 In 2 » -2.07944 
2 x/ 2r (3)4 



(5.25) 



7T 



0.3231(1) 



18 We reduce the accuracy compare to the bosonic mass deformation case in order to speed-up 
computations. For a selected set of 41 sets (5.23) we improve accuracy by roughly a factor of 10 2 . 
More accurate data does not provide noticeable quantitative difference. 
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We do not find a signature for the existence of p cr i t — it appears for supersymmetric 
deformation the boundary value problem (5.20) always have a solution 19 . 

Fig. 5 represent a h (^) (the left plot) and a numerical verification of the first law of 
thermodynamics (3.35) (the right plot): 



<>! H2(P,tt)-f^l +(-6a.tt)-«* 1 ° K) 



green k \ ^£ / J blue I V ^£ / J black 

(da h (£) \ 



(5.26) 



+ 



red 



where subscripts refer to the color of the points on the right plot in Fig. 5. We find 
that (5.26) holds with an accuracy of 1CT 3 . 

5.3 Generic deformation 

Though in this paper we study two special cases of the thermal PW flows — one corre- 
sponding to a bosonic mass deformation where the supergravity scalar x( x ) = 0> an d 
the other one corresponding to a supersymmetric mass deformation where the ratio of 
the coefficients of the non-normalizable modes of the supergravity scalars {p(x), x( x )} 
is the same as for the supersymmetric PW flow — the numerical methods developed 
here can be applied for a generic mass deformation as well. In fact, the Sturm-Liouville 
problem one has to solve in a generic case is exactly the same as in the supersymmetric 



19 Recall for a bosonic deformation p crit w —0.0358(2) which is well inside the set of £ considered in 
Fig. 4. 
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Figure 5: The left plot represents for the supersymmetric mass deformation. 

The red/blue/black points on the right plot represent a' h (t;), 2(p 10 — £p'io(0) an d 
— (Xio(0 _ £%'io(£)) correspondingly. The green points on the right plot represent 
numerical verification of the first law of the thermodynamics, see (3.35). 

case where the boundary conditions are given by (5.20). The only difference from the 
supersymmetric mass deformation case is the (fixed) relation between the coefficients 
of the non-normalizable modes of the supergravity scalars {p(x), x( x )} given by (3.14) 

Recall (see Fig.l) that when m/ = there is a critical value p cr it, such that for 
£ > Peru there is no nonsingular thermal PW flow. On the other hand (see Fig.4), we 
do not find the evidence for p crit for the supersymmetric mass deformation rrif = m^. 
Thus we expect that 



for rrif = implies destabilization of the thermal plasma at low temperatures. What 
we appear to be finding is that this particular instability is cured by turning on rrif ^ 0. 
Verifying the latter conjecture is beyond the scope of this paper. However, we would like 
to point out a different type of instability in a system closely related to the one studies 
here which is eliminated by turning on sufficiently large rrif. Consider thermal PW flow 
with the flat R 3 in (2.7) compactified on a three torus 20 T 3 . This represents holographic 

20 Since a three torus is locally the same as R 3 , such a compactification does not change any of the 
above computations. 




(5.27) 



such that p cr u — > oo as 



A* < 1. As we show in subsection 6.1 the existence of p cr u 
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Figure 6: The left plot represents ratio of thermal PW geometry horizon curvature to 
its curvature at the boundary for bosonic mass deformation. The right plot represents 
this ratio for the supersymmetric mass deformation. 

dual to thermal M = 2* plasma on finite volume T 3 . Analysis identical to that in 
[38] shows that this system is generically non-perturbatively unstable with respect 
to a spontaneous creation of the D3D3 brane pairs 21 . This instability is suppressed 
provided 22 

2 

TYl 

xl > 6pn > ?? 2 >6£>0 A = -4 > 1 . (5.28) 

m b 

So, for the nonperturbative instability due to D3D3 spontaneous brane pair creation 
for the thermal PW flow compactified on T 3 we find that A* = 1. 

5.4 Validity of the supergravity approximation 

Our analysis of the holographic dual to thermal M = 2* gauge theory plasmas are 
done in the supergravity approximation to type lib string theory. In the planar limit 
considered here the string loop corrections are suppressed and the supergravity ap- 
proximation is valid provided curvature invariants of the thermal PW flow geometry 
(2.7) are small. We find that along the flow from the boundary to the horizon the 
absolute value of the Ricci scalar of the five dimensional thermal PW geometry (2.7) 
increases. As a quantitative criteria for such an increase we propose the ratio of the 

21 Similar instability was first identified in [39]. 

22 See eq.(4.16) of [38] where one has to remove (—2) from the inequality (coming from the curvature 
of the 3-space) and replace pn of [38] with 4pn due to a different choice of a radial coordinate. 
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scalar curvature at the horizon to the scalar curvature near the boundary 



k = 




horizon 



(5.29) 



where we used (2.5) and the fact that V boundary = -§. Values 23 of re(f ) for both bosonic 
and supersymmetric mass deformations are presented in Fig. 6. Notice that for the 
'blue' phase 



and thus supergravity approximation for the blue phase of the bosonic mass deforma- 
tion for small values of £ is not valid. 

6 J\f = 2* plasma at finite temperature 

In this section we discuss thermodynamics of the Af = 2* gauge theory plasma for two 
special cases: 

■ bosonic deformation, when to/ = and mj ^ 0; 

■ supersymmetric deformation, when mj = m b = m. 

6.1 Thermal PW flows with m/ = 0, m b ^ 

As explained in subsection 5.1, for the bosonic deformation we observe two phases 
('red' phase and 'blue' phase) of the thermal gauge theory plasma for mass deformation 
parameter £ e (0, p cr u) characterized by two different values for the bosonic condensate 
Pio(0- These are the red and the blue sets of points in Fig. 1. For large values of the 
bosonic deformation parameter (£ > p crit in Fig. 1) thermal Af = 2* gauge theory 
plasma with bosonic deformation discussed does not exist. The set of black points in 
Fig. 1 corresponds to £ < 0, which following (3.17) corresponds to m 2 < 0. As we 
argued in subsection 5.4, we can trust supergravity approximation for the red set of 
points (everywhere) and for black points over the range we studies (corresponding to 
m 2 /T 2 < —15). The supergravity approximation for the blue set of point is much less 
reliable 24 , especially for small values of £. 

23 In order to evaluate r p honzon besides a/j(£) we retained horizon values of the supergravity scalars 
{p(x),x(x)} while solving the Sturm-Liouville problem (5.4) or (5.20). 

24 The smallest absolute value horizon curvature for the blue points is the same as the largest horizon 
curvature for the red points. 



lim k\ 



bosonic 
blue 



(0 = +oo , 



(5.30) 
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•Fbosonicj (g^" ^ T ) 

0.4[ 




= "'<■ 



3 x c 2 



bosonic 



J 



Figure 7: The left plot represents bosonic mass deformation free energy Tbosonic as a 
function of fi = j^. The critical mass is at jj crit ~ 5.28(5). The right plot represents 
speed of sound squared c 2 s bosonic for bosonic mass deformation as a function of fi. Recall 
that the speed of sound in the conformal plasma c 2 sCFT 



l 

3 • 



We use (3.17), (3.27), (3.29) and (3.30) to evaluate the free energy density and the 
speed of sound for the bosonic mass deformation. For the free energy density we set 
an overall (arbitrary) constant part of the free energy 25 



^bosonic = Q (g-Q 

Notice that while in the supersymmetric case one can determine jFg nsj/ unambiguously 
(see (3.36)) by requiring the jF susy (T = 0) = 0, we can not do so for the bosonic 
deformation: as we describe shortly, the bosonic deformation becomes unstable both 
thermodynamically and dynamically at low temperatures. Results of the analysis are 
presented in Figs. 7-8. 

The left plot in Fig. 7 describes the free energy density for the bosonic deformation 
J-bosonic relative (up to a sign) to the free energy density of the conformal gauge theory 
plasma 



CFT 



— 7T 2 iV 2 T 4 



as a function of 



(6.2) 



(6.3) 



Notice that the blue phase has always a larger (or equal) free energy than the red phase 

•F~bosonic,red^y[L) — ^bosonicfiluei^f^) ; (^"^) 



25 Such ambiguity is present both in thermal quantum field theory and in the framework of the 
holographic renormalization of the string theory duals to thermal gauge theories. 
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4 . 9 




Figure 8: More detailed plot of the speed of sound squared c 2 bosonic for bosonic mass 

deformation as a function of p. The speed of sound vanishes at p* ~ 5.40(9). Notice 

that though p* > p cr u, corresponding to p* value of £* ~ 0.035(1) is less than p crit « 

0.03582(2), see Fig. 1. Recall that the speed of sound in the conformal plasma c 2 CFT = 
i 

3" 



with equality achieved only at p = p crit corresponding to p crit according to (3.17) 
(6.3). For the red phase 

~F'bosonic,red(p J 0) 



§7T 2 iV 2 T 4 



(6.5) 



in agreement with the fact that for m& = (or at very high temperatures T — > oo) the 
free energy density must be that of the conformal gauge theory. 

The right plot in Fig. 7 and a more detailed plot in Fig. 8 present the speed of 
sound squared for the thermal M = 2* gauge theory plasma, relative to the speed of 
sound in the conformal gauge theory plasma 

C l,CFT = \ , (6-6) 

for bosonic deformation as a function of fi. Although bosonic deformation phase dis- 
cussed here does not exist for £ > p crit) bosonic phase actually extends to values of 
Hcrit < < At* (see Fig. 8). The reason for this is that relation between £ and fi as 
given by (3.17) involves a factor e 6ah ^ and is not monotonic around £ ~ Pent- A value 
of £* corresponding to p* is actually less than p crit : 

p* w 5.40(9) > p crit w 5.28(5) C w 0.035(1) < p crit w 0.03582(2) . 

(6.7) 
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We find that the speed of sound for the blue phase is superlunimal, which is yet another 
indication that this phase is unphysical. The speed of sound diverges for the blue phase 
as one approach p cr jt either from above or below. The speed of sound for the red phase 
at p = is that as in the CFT plasma; it vanishes at p — p* and becomes purely 
imaginary for £* < £ < p crit (still for the red phase and well within supergravity 
approximation) 

C s,bosonic,red(f Jj ~ 0) — n i C s,bosonic,red(f J ' — P ) — > 

, 3 N (6-8) 

Clbosomcred (p = ^TrV^) < , f < £ < Pcrit ■ 

We would like to stress that conclusions (6.8) are very robust and in fact, immediately 
follow from the existence of p cr u for the bosonic condensate pio(£)- Indeed, existence 
of Pcrit) z.e., impossibility to determine pio(£) for £ > p cr u, along with the assumption 
for the smoothness of pio(£), implies that Pio(0 is not a single valued function of £ 
around p CT u\ moreover, 

c lim p'i (0 = Too, (6.9) 

with the minus sign in the red phase and the plus sign in the blue phase (precisely 
what is presented in Fig. 1). But then from (3.30) 

lim S feosonic (£) = =foo , (6.10) 

?— > Pcrit — 

again with the minus sign in the red phase and the plus sign in the blue phase. Thus 
using (3.27) we see that in the red phase 

lim (s x cl bosonic (0) = -3 , (6.11) 

(in agreement with Figs. 7-8), and for the blue phase 

lim ( 3 X c\ bosonic (£) ) = ±00 , ^bosonicblue^oo) = 3 , (6.12) 

where < £oo < p cr u- Since in the red phase at £ = the speed of sound squared is | > 
and given (6.11) there is a guaranteed value of £* < p cr u at which the speed of sound 
vanishes, and for £* < £ < p crit the speed of sound is purely imaginary. As emphasized 
in [35], a thermal system with c 2 s < is unstable simultaneously thermodynamically 
(it has a negative specific heat) and dynamically (amplitude of small pressure/energy 
density fluctuations exponentially grows). 
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fsusy — ^Fsusy/ (g 71 " N T ) „ ^(fsusy) 




Figure 9: The left plot represents supersymmetric mass deformation free energy fF susy 
as a function of The right plot represents ln(f susy ) as a function of ^ (black points) 
and a straight line fit (red curve) for the last 2500 points of the distribution. 



Finally, we do not see thermodynamic instability in thermal Af = 2* plasma with 
bosonic mass deformation for m 2 < 0, see Fig. 7. In fact, we find that for negative 
values of \i the speed of sound slowly (and monotonically) decreases — it reaches about 
90% of the conformal plasma speed of sound (6.6) at 26 /i ~ —35. Of course, this does 
not exclude the fact that thermal Af = 2* plasma with m 2 < might have dynamical 
instability without accompanying thermodynamics instability [35]. 

6.2 Thermal PW flows with rrif = = m 

Bosonic and fermionic condensates {pio(0> Xio(0} f° r thermal Af = 2* gauge theory 
plasma with supersymmetric mass deformation are presented in Fig. 4. As we argued 
in subsection 5.4, supergravity approximation here is valid of the whole range of de- 
formation parameter £ considered. We use (3.17), (3.27), (3.33), (3.34) and (3.36) to 
evaluate the free energy density and the speed of sound for the supersymmetric mass 
deformations. Results of the analysis are presented in Figs. 9-10. We do not see any 
signature of the phase transition 27 , or the thermodynamic instability. 

Probably the most unexpected result of the analysis is that strongly coupled ther- 
mal Af = 2* plasma for supersymmetric mass deformation behaves much closer to a 
conformal plasma than one would naively expect. The right plot on Fig. 9 presents 
the logarithm of the ratio of the Af = 2* plasma free energy density to the free energy 

26 This is outside the range of the right plot in Fig. 7. 

27 This would show up as the free energy T susy (^t) changing sign for a particular value of ^. 
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Figure 10: The speed of sound squared c 2 snsj/ for supersymmetric mass deformation as 



a function of ^. 

density of the conformal plasma (6.2) 

Mf susy ) = \n^L, (6.13) 

J~CFT 

as a function of ^. We did a straight line line fit (solid red line on the right plot on 
Fig. 9) for the low temperature 28 



FIT 



(x) =0.118(8) -0.138(1) x. (6.14) 



Thus, to a good approximation the free energy of the thermal Af = 2* plasma can be 
approximated by 

TtZ y ° x ~ 1-1(3) x Tcft e'^Tr oc _ T 4 e -f^ ? ( 615 ) 

so that noticeable deviation from conformality occurs at T ~ |m. 

It is instructing to compare (6.15) with the zero coupling result for the M = 2* 
plasma. At zero coupling we have 29 N 2 massless M = 2 vector multiples and N 2 
massive M = 2 hypermultiplets. The free energy density of the massless states is just 
\Fcft i n (6-2), while the free energy density of a mass m spin s particle is given 
by [40,41] 

r = J»^±?£KJ%), (, 16) 



28 In practice, wc define low temperatures as the subset of the mass deformation parameters ^ >^ f . 

29t 



We neglect (—1) in the large N limit. 
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where r\ — ±1 for bosons/fermions correspondingly. As low temperatures | > 1 we 
can use asymptotic expansion for the Bessel function K 2 (x) 

K ^ x) ^^c e ~ X ' a:>1 ' (6 ' 17) 

to conclude that 

(2s + l)m 3 / 2 T 5 / 2 m m , 

^ ( 2 V e ~"> T >>1 - (6 - 18) 

Altogether we find that the free energy of the thermal Af = 2* plasma for supersym- 
metric mass deformation at zero coupling is 

tt 2 N 2 T 4 ( 24V2 /m\3/2 m \ m 
= — L2- I 1 + !W (t) e ~ W ) ' T > 1 • (6 - 19) 

From (6.19), we find that the free energy density is roughly 60% of the high temperature 
result at T ~ m, while the corresponding number for the strongly coupled plasma (6.14) 
is about 98%. 

The hydrodynamic properties and the jet quenching of strongly coupled nonconfor- 
mal gauge theory plasma are typically parameterized in terms of [25,42] S 

6= Q - cf) . (6.20) 

For the strongly coupled Af = 2* thermal plasma at 30 T ~ 1.5m we find (see Fig. 10) 

<^ =2 Vi. 5 m~2x 10- 3 , (6.21) 

which further emphasizes the striking similarity between strongly coupled Af = 2* 
plasma and conformal Af = 4 gauge theory plasma. 



7 Conclusion and future directions 

In this paper we presented detailed analysis of the thermodynamics of the strongly 
coupled Af = 2* gauge theory plasma. We considered two special cases of the mass 
deformations: bosonic mass deformation (where only the bosonic components of the 

30 The temperature of the quark-gluon plasma at RHIC is of order the QCD deconfinement temper- 
ature [43]. 
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Figure 11: The left plot represents bosonic mass deformation energy density Sb osonic as 
a function of ^ . The value ^/Jjl* (see Fig. 8) denotes the ratio ^ for which the speed 
of sound vanishes. The right plot represents supersymmetric mass deformation energy 
density S susy as a function of ^. 

Af = 2 hypermultiplet get nonzero mass) and supersymmetric mass deformation (with 
bosonic and fermionic components of the hypermultiplet getting the same mass). We 
argued that both our theoretical framework and the numerical procedures are under 
control, given a highly nontrivial consistency check of the first law of thermodynamics 
for the extracted data. We also argued that the supergravity approximation is under 
control, as the curvature invariants of the dual geometry remain small. 

The most important result of the analysis is the striking similarity between ther- 
modynamics of the mass deformed Af = 2* gauge theory plasma at temperatures of 
order the mass scale and the thermodynamics of the conformal Af = 4 SYM plasma. 
We would like to reemphasize the latter fact by presenting the plots for the energy 
energy for various mass deformations, see Fig. 11. 

We could like to concluded with some open problems. 

■ In is important to generalize our analysis to other examples of the strongly coupled 
nonconformal gauge theory plasma. Most notably the cascading gauge theory plasma 
[18-20,22]. 

■ It is interesting to evaluate hydrodynamic properties and jet quenching in Af = 2* 
gauge theory plasma in the low temperature regime. 

■ The speed of sound waves in Af = 2* gauge theory plasma for the bosonic mass 
deformation was found to vanish for certain value of ^r. It is interesting to study 
hydrodynamics of the gauge theory plasma in the regime of the vanishing speed of 
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sound and determine the critical exponent of the bulk viscosity. 
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